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ABSTRACT 

We report on the results of the extensive multi-wavelength campaign from optical to GeV 
y-rays of the 2014 periastron passage of PSR B1259-63, which is a unique high-mass y-ray 
emitting binary system with a young pulsar companion. Observations demonstrate the stable 
nature of the post-periastron GeV flare and prove the coincidence of the flare with the start of 
rapid decay of the Hu equivalent width, usually interpreted as a disruption of the Be stellar 
disk. Intensive X-ray observations reveal changes in the X-ray spectral behaviour happening at 
the moment of the GeV flare. We demonstrate that these changes can be naturally explained as 
a result of synchrotron cooling of monoenergetic relativistic electrons injected into the system 
during the GeV flare. 

Key words: gamma rays: stars - pulsars: individual: PSR B1259-63 - stars: emission-line. 
Be - X-rays: binaries - X-rays: individual: PSR B1259-63 


1 INTRODUCTION 

In the binary system PSR B1259-63 a 47.76 ms radio pul¬ 
sar is in a highly eccentric orbit (e ~ 0 .87, P «_3.4 years) 

around the massive 09.5Ve star LS 2883 l ljohnston et alj|l992l : 
iNegueruela et alj 1201 lh . The optical spectrum of the companion 
shows evidence of an equatorial disk (thus this star is gener¬ 
ally classified as a Be star), whic h is th o ught to b e inclin ed w ith 
respect to the orbital plane llMelatos. Johnston & Melroselll995l) . 
The pulsar crosses the disk plane twice each orbit, just before 


and after the periastron passage. The minimum approach between 
the pulsar a n d massive star is about ~ 0.9 astronomical unit 
jNegueruela_etalj 201 lh. whic h is roughly the size of the equato¬ 
rial disk lljohnston et alJlT992h . Interaction between the relativis¬ 
tic pulsar wind and the wind and photon field of the Be star is 
believed to give rise to the observed unpulsed e mission . X- ray 
emission is observed througho ut the orbit (e.g. iHiravama et alj 
Il999l : IChemvakova et alj|2006l) and the unpulsed radio, GeV and 
TeV emission is observed within a few months of periastron pas- 
sage jjohnston et alj 1 19991120051 : iKirk, Ball & Skiaeraasenl 1 19991 : 
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Abdo et al]|201ll : IChemvakova et al1l2014l : lAharonian et al.ll2005l 
20091) . 


The presence of a nearby X-ray pulsar, 
2RXP J130159.6-635806, located only 10 arcminutes away 
from PSR B1259-63 makes it difficult to observe the source with 
a non-imaging instrument, like Suzaku HXD-PIN detectors, or 
instruments with a low angular resolution, like INTEGRAL. For 
these instruments both sources influence the observed emission, 
and in order to reconstruct the spectral shape of PSR B1259-63 it 
is necessary to know the spectrum of 2RXP J130159.6-635806. 
Recent NuSTAR observations of the 2RXP J130159.6-635806 
dKrivonos et alj l2015h allowed us to reconstruct the broad band, 
0.5 - 60 keV X-ray spectrum of PSR B1259-63 as observed by 
Suzaku in 2011, and by INTEGRAL in 2015. 

The most puzzling feature of the orbital modulation of the 
source is the huge flare in the GeV band starting ten days af¬ 
ter the post-periastron passage of the pulsar through the stel- 
lar disk with no obvious counterpart at any other wavelength 
dAbdo et al .120 111 : IChemvakova et alj2014l) . T he origin of this flar e 
was widely discussed i n the literature (see e.g.lAbdo et aid d201 lh: 
IPetri & D ubusl 12011): Khangulvan et aid d20 1 2 I) : Mochol & Kirkl 
d20 1 3l) : lDubus & Ceruttil d20 1 3ll ). but the lack of observational data 
prevented any firm conclusions being made. It was even not clear 
in advance whether the GeV flare would repeat in 2014. 

In this paper we present the results of the extensive multi¬ 
wavelength campaign of the latest PSR B1259-63 periastron pas¬ 
sage, which happened on 2014 May 4 (t J ,=MJD56781.42). One of 
the main aims of this campaign was multiwavelength monitoring 
of PSR B1259-63 around the periastron and the moment of GeV 
flare in order to reveal the true nature of the flare. The paper is or¬ 
ganised the following way. In section[2]we present results of optical 
monitoring with the Southern African Large Telescope (SALT) and 
SAA0 1.9-m telescopes. In section[3]details of X-ray data analysis 
along with the results of X-ray monitoring with SWIFT, NuSTAR, 
and INTEGRAL are given. In this section we also present for the 
first time the broadband X-ray spectrum of PSR B1259-63 as seen 
in 2010 by Suzaku. Section[4]explains details of Fermi-L data anal¬ 
ysis. Finally all the obtained results are discussed in section[5] 


2 OPTICAL SPECTROSCOPY 
2.1 SAA0 1.9-m telescope 

Spectroscopic observations of LS 2883, the primary of the system, 
were performed with the SAA0 1.9-m telescope between the 2014 
lune 6 and 2014 lune 11. The Cassegrain spectrograph was set with 
grating number 5 to achieve a wavelength coverage of 6350-7050 
A with a blaze of 6800 A and a resolving power R = A/AA = 6563 
in the vicinity of the H a line. 

LS 2883 was observed for a total of 3 nights using 30s ex¬ 
posure times. 240 individual spectra were taken on the first night, 
245 on the second and 196 on the third. Copper-Neon comparison 
lamp spectra were obtained throughout each night for wavelength 
calibration. 

The spectra were flat fielded, extracted and wavelength cal¬ 
ibrated using the standard slit spectroscopy routines in IRAF. To 
improve the signal-to-noise (S/N) of weak lines the spectra for each 
night were coadded. The mean spectra were then normalised to a 
unit continuum. 

The equivalent widths of the H a and He 1(6678 A) lines were 
measured by integrating over the emission line profiles. Due to the 


Table 1. Journal of SAA0 1.9-m telescope 2014 observations of the source 
PSR B1259-63 around the time of GeV flare. 


Date 

MJD 

(days) 

t-tp 

(days) 

V 

-w Ha 

(A) 

^IT 66 78 

(A) 

2014-06-06 

56814.74 

33.3 

294.8 

70.28 ± 0.73 

0.44 ± 0.05 

2014-06-07 

56815.77 

34.3 

295.8 

71.48 ±0.81 

0.43 ± 0.05 

2014-06-10 

56818.71 

37.3 

298.4 

69.61 ± 1.16 

0.44 ± 0.05 

2014-06-11 

56819.87 

38.4 

299.3 

71.17 ±4.39 

- 

2014-06-12 

56820.72 

39.3 

300.0 

67.82 ± 4.22 

- 

2014-06-13 

56821.69 

40.3 

300.7 

67.68 ±4.18 

- 

2014-06-16 

56824.72 

43.3 

302.9 

66.39 ±4.16 

- 

2014-06-17 

56825.74 

44.4 

303.6 

65.84 ±4.10 

- 


high S/N of the co-added spectra, the errors in W(Hff) and W(He 
1(6678 A)) are due principally to the placement of the continuum. 
The larger percentage errors in W(He 1(6678 A)) is a result of an 
absorption feature in the vicinity, and the lower line strength. By 
choosing appropriate upper and lower bounds on the continuum 
levels, the error margins on the equivalent widths were calculated. 

Further observations of LS 2883 were made between the UT 
dates 2014 June 11 and 2014 June 17. The Cassegrain spectrograph 
was set with grating number 6 to achieve a wavelength coverage of 
4100-5400 A with a blaze of 4600 A and a resolving power R = 
A/AA = 2430 in the vicinity of the line. 

In this mode LS 2883 was observed for 5 more nights with ei¬ 
ther 1 or 2 spectra taken each night. The exposure times each night 
varied from 350s to 700s depending on the number of spectra taken. 
Copper-argon comparison lamp spectra were obtained before and 
after each sequence for wavelength calibration. 

An attempt was made to infer the equivalent width of the 
line using the equivalent width of the line. Simultaneous 
measurements of the and lines have been previously made 
and the ratio of W(H„) to W(H^) was found to be 11.74 ± 0.67 
dNegueruela et al.ll201 lh . The application of this result to our H/j 
data allowed for the deduction of W(H„) from the measured WJH^) 
values. The larger error margins on the converted H ff data relative 
to the Hj, data from the first three nights of observations is due 
to a combination of larger percentage errors on the WJH^) values 
and the propagat ion of the uncertainty in the W(H„)/W(H/j) relation 
from lNegueruela et al. (201jj). _ 

The data of Negueruela et al] d20 1 lh were taken at a phase of 
0.62, while our observations are made at a phase of 0.02. The as¬ 
sumption that the W(H a )/W(Hyj) ratio remains constant with phase 
is not completely justified considering that the physical conditions 
of the disk drastically varies as the pulsar moves through. However, 
spectroscopic observations from SALT were made simultaneously 
for two of the nights. The agreement between the two data sets sug¬ 
gests that the W(H <r )/W(H j g) relation does not change with phase 
by more than the meas ured uncertainty in the equivalent widths of 
Negueruela et al] J201 ll) . 


2.2 SALT observations 

As discussed in van Soelen et al. (in prep) optical spectroscopy was 
undertaken with the South Afric an Large Telescop e (SALT) using 
the Robert Stobie Spectrograph dBurgh et ai1l2003l) . between 2014 
April 30 and 2014 July 21. The RSS was used in a long slit mode 
with a wavelength coverage of 6176.6 - 6983.0 A, with a resolu¬ 
tion of R = 11021 at the central wavelength. A typical observation 
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Figure 1. Orbital light curves of PSR B1259-63 around periastron for several passages. Dashed lines correspo nd to the periastron and to the moments of 
disappearance (first non-detection) and reappearance (first detection) of the pulsed emission, as observed in 2010 lAbdo et alj201 ll) . Dotted line corresponds 
to the first appearance of the detection in GeV band at a day time scale. Panel a: Fermi-L. flux measurements in the E > 100 MeV energy range for the 2010 
and 2014 periastron passages with a weekly bi n size . Flux is given in 10 6 cm~ 2 s 1 . Panel b: X-ray fluxes from 2014 (this work) and three last periastron 
passages l.Abdo etaOl201 ltlChemvakova et al .11 20091) . Flux is given in 10 11 erg cm -2 s 1 . The typical error of the X-ray data is smaller than the size of the 
symbols. Panel c: Evolution of the equivalent widths of Her. 


consisted of 3-4 camera exposures (total of ~ 476 to 500 s) which 
were co-added to achieve a higher SNR. 

Data were reduced following the standard IRAF procedures, 
and flux shape correction was done using the spectroscopic stan¬ 
dard LTT4364 (observed on 2014 May 11). 

The equivalent width of the Her line was measured by integrat¬ 
ing over the line profile (within IRAF) with an assumed linear con¬ 
tinuum between the selected positions. The statistical error in the 
equivalent width has bee n esti mated using the method discussed in 
IVollmann & Eversberg ( 2006 ). In order to include an estimate of 
the uncertainty due to the line continuum, each measurement was 
performed more than once and any variation in the answer was in¬ 
cluded in the error estimate. In general this was smaller than the 


statistical estimate, though a more rigorous estimate of the error 
introduced by continuum placement (such as used for the SAAO 
1.9-m) will slightly increase the error. 


Since there was limited overlap between the SALT and SAAO 
1.9-m observations (and this is for the SALT data which suffered 
from instrumental problems), we performed a comparison between 
the analysis methods used for the different data sets by analysing a 
selection of spectra from both the SALT and SAAO 1.9m observa¬ 
tions using both analysis methods, and confirmed that both methods 
gave results that agreed within the uncertainty. The results are listed 
in Table [2] 
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Table 2. Journal of SALT telescope 2014 observations. 


Date 

MJD 

(days) 

t tp 

(days) 

V 

-w Ha 

(A) 

2014-04-01 

56748.0 

-33.4 

65.06 

54.88±0.68 

2014-04-05 

56752.9 

-28.5 

70.38 

55.27±0.72 

2014-04-11 

56758.0 

-23.4 

77.48 

55.56±0.62 

2014-04-16 

56763.9 

-17.5 

88.89 

58.59±0.10 

2014-04-25 

56772.9 

-8.5 

119.94 

55.37±0.55 

2014-04-29 

56776.8 

-4.6 

143.18 

56.78± 0.69 

2014-04-30 

56777.9 

-3.5 

151.17 

56.55± 0.59 

2014-05-01 

56778.8 

-2.6 

158.13 

57.97± 0.54 

2014-05-03 

56780.8 

-0.7 

174.59 

60.90±1.02 

2014-05-11 

56788.8 

7.4 

233.86 

68.46± 0.55 

2014-05-17 

56794.8 

13.4 

259.57 

72.70± 0.71 

2014-05-23 

56800.8 

19.4 

275.12 

71.32± 0.71 

2014-05-31 

56808.8 

27.4 

288.14 

70.62±2.50 

2014-06-06 

56814.7 

33.3 

294.77 

68.59 ±4.39 

2014-06-12 

56820.7 

39.3 

299.99 

64.96± 0.70 

2014-06-13 

56821.7 

40.3 

300.75 

63.70± 0.54 

2014-06-16 

56824.7 

43.3 

302.88 

60.83± 0.51 

2014-06-17 

56825.8 

44.4 

303.61 

61.46± 0.75 

2014-06-19 

56827.7 

46.3 

304.82 

60.54± 0.70 

2014-06-21 

56829.7 

48.3 

306.01 

60.12± 0.82 

2014-06-22 

56830.8 

49.4 

306.64 

60.14± 0.65 

2014-06-29 

56837.8 

56.3 

310.21 

60.15±3.21 

2014-07-09 

56847.8 

66.4 

314.37 

57.39±0.79 

2014-07-13 

56851.7 

70.3 

315.77 

55.89±0.49 

2014-07-21 

56859.7 

78.3 

318.33 

56.69± 0.55 


3 X-RAY OBSERVATIONS AND RESULTS 

We conducted an X-ray monitoring campaign on PSR B1259-63 
with Swift , NuSTAR and INTEGRAL telescopes, covering the pe¬ 
riod between f p - 15 days and f p + 65 days. These observations 
are summarized in Tables as and [6] Tables [3] and [4] list iden¬ 
tifiers for the data set, UT date, MJD, time relative to periastron 
passage and exposure time for each observation. In this paper we 
also use new NuSTAR observations of the nearby X-ray pulsar 
2RXP J130159.6-635806 to reconstruct the broad 0.5 - 50 keV 
X-ray spectrum as observed by Suzaku in 2014. All 2014 observa¬ 
tions can be well fitted with an absorbed power law model. 


3.1 Swift/XRT data reduction 

The Swift observatory jGehrels et al.ll2004 ) provides the possibil¬ 
ity to monitor sources of X-ray emission on very different time 
scales. In this work we use observations covering almost three 
months around periastron passage of PSR B1259-63 between 2014 
April 20 and July 8. In Table [3] a log of the used observations is 
given. The data were processed using tools and packages available 
in FT00LS/HEAS0FT 6.14. 

XRT observed PSR B1259-63 both in Photon Counting (PC) 
and Windowed Timing (WT) modes. Initial cleaning of events has 
been done using xrtpipeline with st andard parameters. The fur¬ 
ther analysis was performed following lEvans et al j 12009 ). In par¬ 
ticular, in the PC mode the source extraction region was a cir¬ 
cle with radii from 5 to 30 pixels depending on the count rate 
fevans et al .l2009h : in the WT mode radius of the source extraction 
region was 25 pixels. The background was collected over the an¬ 
nulus region with an inner (outer) radius of 60 (110) pixels in both 
observational modes. The count rate from the source was too low to 
pile up the detector in all observation except one (00030966025), 


where we excluded the inner region of the source aperture with a 
radius of 4 pixels. 

The obtained spectra were grouped to have at least 1 count 
bin -1 using the FTOOLS grppha. To avoid any problems caused 
by the calibration uncertainties at low energiesU, we restricted our 
spectral analysis to 1.0 - 10 keV. The errors reported in this work 
are purely statistical and correspond to a lcr confidence level. 

To estimate the 68% confidence ranges for the unabsorbed 
fluxes, we computed the C-statistic likelihood profile on the grid 
of column density ( n H ) and power law slope (T) values. According 
to the Wilks theorem, the likelihood values are distributed around 
the true minimum with a^ 2 distribution with two degrees of free¬ 
dom, since we have two free parameters. We thus converted the 
obtained likelihood values to probabilities and performed Monte 
Carlo simulations in order to convert the computed likelihood pro¬ 
file to the distribution of the unabsorbed fluxes. For each of those 
distributions (for every single observation) we calculated the mean 
and its (asymmetric) confidence range by computing the limiting 
flux values below and above the mean, that contain between them 
68% of all the simulated fluxes. 


3.2 NuSTAR observations and data analysis 

Previous X-ray observations of PSR B1259-63 suffered due to 
a lack of high quality X-ray data above 10 keV. The first hard 
X-ray measurement of the system was done with the OSSE in¬ 
strument in the 40 - 500 keV band durin g monitoring of the 
1994 periastron passage l lGrove et' aD 1 19951) . The OSSE experi¬ 
ment could not provide spatially-resolved flux measurements of 
the source due to its non-imaging design and, therefore, flux pollu¬ 
tion from nearby sources could not be fully excluded. The first at¬ 
tempt to perform hard X-ray imaging of PSR B1259- 63 was done 
with INTEGRAL during periastron passage in 2004 iShaw et aD 
120041) . which demonstrated the importance of spatially-resolved 
observations. It was shown that the variable hard X-ray source 
2RXP J130159.6-635806 located 10' away from PSR B1259-63, 
significantly contributes to the hard X-ray flux measured from this 
region. Despite the progress made in understanding the contami¬ 
nation from 2RXP J130159.6-635806 INTEGRAL was unable to 
provide non-contaminated X-ray measurements of PSR B1259-63 
due to its 12' angular resolution. Later periastron observations with 
Suzaku were also not able to spatially separate hard X-ray emis¬ 
sion from PSR B1259-63 and 2RXP J130159.6-635806, forc¬ 
ing authors to assume their hard X-ray spectral shapes based on 
2-10 keV data JUchivama et al. I [20091) . Suzaku observations in¬ 
dicated the possible break during the first interaction of the pulsar 
with the disk. 

Significant progress has been recently made thanks to the 
launch of the NuSTAR mission (Harrison et al., 2013), the first 
hard X-ray imaging telescope in space working at energies above 
10 keV. NuSTAR team initiated an observational campaign during 
periastron passage of PSR B1259-63 in 2014. In total NuSTAR per¬ 
formed five observations for a total of 146 ks, targeting the source 
on-axis of the two co-aligned X-ray telescopes. The 10 m focusing 
mirror and detector modules provide 58' (half-power diameter) and 
18' (FWHM) imaging resolution over the 3-79 keV X-ray band, 
with a spectral resolution of 400 eV (FWHM) at 10 keV. The NuS¬ 
TAR field of view (FOV) is 12' x 12' at 10 keV as defined by the 


1 http://www.swift.ac.uk/analysis/xrt/digest_cal.php 
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Table 3. S wift/XRT observations of the source PSR B1259-63 


Data 

Set 

Obs Id 

Date 

Date, 

MJD (UTC) 

t ~ tp, 

d 

Exposure, 

s 

XRT count rate' 7 , 
cts s -1 

XRT 

mode 

Swl3 

00030966018 

2014-04-20 

56767.1736 

-14.561 

1171.2 

0.36 

PC 

Swl4 

00080099001 

2014-04-20 

56767.4857 

-14.249 

1740.6 

0.44 

PC 

Swl5 

00030966019 

2014-04-27 

56774.9560 

-6.779 

4015.7 

0.29 

PC 

Swl6 

00030966020 

2014-05-04 

56781.9594 

0.224 

4075.6 

0.22 

PC 

Swl7 

00030966021 

2014-05-06 

56783.8933 

2.158 

1470.9 

0.20 

PC 

SwlS 

00030966022 

2014-05-07 

56784.0872 

2.352 

2404.9 

0.19 

PC 

Swl9 

00030966023 

2014-05-09 

56786.1540 

4.419 

1475.9 

0.16 

PC 

Sw20 

00030966024 

2014-05-17 

56794.5981 

12.863 

179.8 

0.44 

PC 

Sw21 

00030966025 

2014-05-19 

56796.8441 

15.109 

1987.8 

0.61 

PC 

Sw22 

00030966026 

2014-05-22 

56799.4340 

17.699 

2877.3 

1.00 

WT 

Sw23 

00030966027 

2014-05-25 

56802.4956 

20.761 

2925.1 

0.74 

WT 

Sw24 

00030966028 

2014-05-28 

56805.6925 

23.958 

2991.7 

0.65 

WT 

Sw25 

00030966029 

2014-06-02 

56810.9245 

29.189 

1504.4 

0.51 

WT 

Sw26 

00030966030 

2014-06-03 

56811.0620 

29.327 

970.0 

0.12 

WT 

Sw27 

00030966032 

2014-06-07 

56815.0953 

33.360 

2696.4 

0.44 

WT 

Sw28 

00030966033 

2014-06-11 

56819.2927 

37.558 

3983.2 

0.54 

WT 

Sw29 

00030966035 

2014-06-12 

56820.4606 

38.726 

999.6 

0.54 

WT 

Sw30 

00030966036 

2014-06-13 

56821.8494 

40.114 

2036.3 

0.51 

WT 

Sw31 

00030966037 

2014-06-14 

56822.8488 

41.114 

1834.4 

0.49 

WT 

Sw32 

00030966038 

2014-06-15 

56823.2857 

41.551 

767.3 

0.44 

WT 

Sw33 

00030966039 

2014-06-16 

56824.5833 

42.848 

2026.2 

0.48 

WT 

Sw34 

00030966041 

2014-06-18 

56826.2590 

44.524 

981.5 

0.50 

WT 

Sw35 

00030966043 

2014-06-25 

56833.3399 

51.605 

3284.0 

0.28 

PC 

Sw36 

00030966044 

2014-06-26 

56834.5113 

52.776 

1513.4 

0.27 

PC 

Sw37 

00030966047 

2014-06-27 

56835.4127 

53.678 

2559.7 

0.25 

PC 

Sw38 

00030966048 

2014-07-02 

56840.4456 

58.711 

2964.3 

0.26 

PC 

Sw39 

00030966049 

2014-07-06 

56844.7580 

63.023 

3264.0 

0.20 

PC 

Sw40 

00030966050 

2014-07-07 

56845.2536 

63.519 

3785.9 

0.17 

PC 

Sw41 

00030966051 

2014-07-08 

56846.5879 

64.853 

3965.7 

0.18 

PC 


° Total count rate in 0.5-10 keV energy range 


full width at half intensity. The nominal reconstructed coordinates 
are accurate to 8" (90% confidence level) (Harrison et al. 2013). 

We reduced and analyzed data using the NuSTAR Data Anal¬ 
ysis Software (NuSTARDAS) vl.3.1, which is part of the HEA- 
SOFT 6.15 package. The data were filtered for intervals of high 
background and corrected for angular offset using catalogued coor¬ 
dinates of PSR B1259-63. 

During most of the NuSTAR observations 
2RXP J130159.6-635806 was out of the FOV, and only in 
one of them (ObsID: 30002017004) it appeared in the corner at 
a large off-axis distance, which is demonstrated in Fig. [2] It is 
clearly seen from the Fig. [2] that NuSTAR fully resolved flux from 
PSR B1259-63 and 2RXP J130159.6-635806 allowing one to 
conduct for the first time hard X-ray observation of PSR B1259-63 
around periastron without any contamination from the nearby 
variable and hard 2RXP J130159.6-635806 source. As discov¬ 
ered by IChernvakova et aD ( l2005h . 2RXP J130159.6-635806 
shows coherent spun-up pulsations at a period of ~ 700 s. It 
is interesting to note that serendipitous NuSTAR detection of 
2RXP J130159.6-635806 at the comer of the FOV allowe d 
us to measure its flux pulsations period iKrivonos et alj 120151) . 
which is in full agreemen t with spin-up rate evolution reported in 
IChernvakova et all l l2005t) . 

We performed NuSTAR spectral analysis using nuproducts 
vO.2.5, which is part of NuSTARDAS, to generate NuSTAR re¬ 
sponse matrix (RMF) and effective area (ARF) files for an on-axis 
point source. To extract the source spectrum we utilized circular 
region centred at the source position with 70" radius which com- 


Unfolded Spectrum 



5 10 20 50 

Energy (keV) 


Figure 3. Spectra of NuSTAR NST2 (black and red) and NST3 (green and 
blue) observations. 


prises ~ 80% of the encircled energy (Harrison et al., 2013). The 
background spectrum was extracted from a 115" circular region 
positioned away from the source, on the same chip (each detec¬ 
tor module has four chips) avoiding chip edges and gaps between 
them. 
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Table 4. NuSTAR observations of the source PSRB1259-63 


Data 

Set 

Obs Id 

Date 

Date, 

MJD (UTC) 

t-t p , 

d 

Exposure, 

s 

NuSTAR count rate", 
(FPMA. FPMB) 10 1 cts s“‘ 

NSTl 

30002017002 

2014-04-20T05:06:07 - 2014-04-20T22:l 1:07 

56767.23 

-14.189 

30517.85 

6.00,5.51 

NST2 

30002017004 

2014-05-04T10:01:07 - 2014-05-05T04:36:07 

56781.42 

0.005 

33317.07 

3.25,2.92 

NST3 

30002017006 

2014-05-28T10:11:07 - 2014-05-28T22:31:07 

56805.44 

24.021 

27089.72 

7.15,6.69 

NST4 

30002017008 

2014-06-02T 19:21:07 - 2014-06-03T08:06:07 

56810.82 

29.405 

26288.55 

6.15,5.60 

NST5 

30002017010 

2014-06-14T 17:21:07 - 2014-06-15T10:06:07 

56822.73 

41.312 

29106.89 

5.56,5.25 


a Total count rate in 3 - 79 keV energy range 



Figure 2. NuSTAR 3-79 keV image of the region around PSR B1259—63 (in the centre) obtained with two focal plane modules FPMA (left) and FPMB 
(right). This shows the one out of five observations performed by NuSTAR where 2RXP J130159.6—635806 is visible in the lower right corner. The images are 
smoothed with 7" Gaussian kernel for convenience. 


3.3 Reanalysis of 2010 Suzaku observations. 


PSR B1259-63 was extensi vely monit ored by Suzaku during it s 
2007 l lUchivama et al.lEoQ9h and 2010 JChernvakova et aD £2014 ). 
observations Sz9 - Szll in Table [5} periastron passages. Suzaku 
2011 broad spectrum of PSR B1259-63 was never presented in 
the literature. The Suzaku observation s we re perfo rmed with the 
X-ray Imaging Spectrometer (XIS: Ifevama_et_al] j2007f)) in 0.3 - 
12 keV and the hard X-ray detector (HXD; Takahashi et alJd2007l) ) 
in 13 - 600 keV. The FOV of both XIS and HXD also contains 
the X-ray source, 2RXP J130159.6-635806. This X-ray pulsar has 
a luminosity comparable to the one of PSR B1259-63 during the 
periastron passage and makes significant contribution to the hard 
X-ray flux measured by the non-imaging HXD-PIN detector. Re¬ 
cent NuSTAR observations of 2RXP J130159.6-6 35806 allowed 
for the first time to measure its broadband spectrum l lKrivonos et alj 
120151 ). The spectrum of 2RXP 1130159.6-635806 is well described 
by the absorbed power law with a high energy cut off, F(e ) = 
Xe~ r exp(-(e - e c )/e f ). 


To reconstruct the 2011 broadband spectrum of PSR B1259- 
63 we reanalyzed historical Suzaku data using the HEASOFT soft¬ 
ware package (version 6.16), with calibration files distributed on 


01.07.2014. For the XIS and HXD-PIN, we made use of cleaned 
event files, in which standard screening was applied. 

HXD data is a sum of emission coming from PSR B1259-63 
and 2RXP 1130159.6-635806. In our analysis we assumed that the 
shape of the 2RXP J130159.6-635806 spectrum was the same in 
2010 and 2014, and only normalization changed from one observa¬ 
tion to another. The normalization of the 2RXP J130159.6-635806 
spectrum was defined from the simultaneous observations of the 
imaging instruments, PIN/Suzaku for SzlO, and MOSl,2 /XMM- 
Newton for Sz9 and Szll. The s e XMM -Newton observations are 
described in I Chernyakova et alj j2014l) . Simultaneous fits of Sz9, 
SzlO, Szll and NuSTAR data give a good fit (y 2 = 6795.64 for 
6686 degrees of freedom) with the following parameters: Nh = 
2.73±0.07, e c = 6.05 ± 0.20, e f = 13.94 ± 0.46, T = 1.14 ± 0.03, 
see Figure[4] The resulting parameters for PSR B1259-63 and 2 - 
10 keV absorbed fluxes for 2RXP 1130159.6-635806 are given in 
Table |5] 

3.4 INTEGRAL Observations and data analysis 

INTEGRAL IWinkler et alj|2003l ) is a y-ray mission covering the 
energy the range 15 keV - 10 MeV. Observations are carried 
out in individual Science Windows (ScW), which have a typi- 
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Table 5. Journal of 2011 Suzaku observations of PSR B1259-63. 


Data 

Set 

Date 

MJD 

(days) 

t tp 

(days) 

<t> 

(deg) 

Exposure 

(ks) 

N h 

10 22 cm' 2 

r 

Fps*5(l-10keV) 

10“ 11 erg cm -2 s -1 

F^(2-10keV) 

10“ 11 erg cm -2 s -1 

Sz9 

2011-01-05 

55566.8 

22 

99.6 

90.0 

0.54±0.01 

1.78±0.01 

2.84±0.01 

2.10±0.03 

SzlO 

2011-01-24 

55585.6 

41 

121.0 

40.3 

0.49±0.01 

1.54±0.01 

1 .68±0.02 

2.05±0.03 

Szl 1 

2011 -02-02 

55594.2 

49 

126.5 

21.5 

0.46±0.02 

1.46±0.02 

1.47±0.02 

2.26±0.03 




Figure 4. Suzaku XIS (blue crosses) and PIN spectra (green crosses) for S9 (left), S10 (middle) and S1I (rights) observations of 2010. The PIN spectrum 
after background subtraction is decomposed into the two components through joint fitting: PSR B1259-63 power law model (dashed line), and 2RXP J1301 
cutoff power law model (dash-dotted line). 2RXP J130159.6—635806 was in the field of view only during S10 observation (red crosses). For the reconstruction 
ofPSR B1259-63 PIN spectrum for S9 and SI 1 observations we have used simultaneous XMM-Newton data (red crosses on the left and right pictures). 


cal time duration of about 2000s. In 2014 INTEGRAL observed 
PSR B1259-63 periastron passage from June 26 to July 10. The 
data cover revolutions 1429-1433, adding up to a total effective 
exposure time of 141 ks for IBIS/ISGRI in 18-60 keV and 259 
ks for JEMX1 & JEMX2 combined. The data were analyzed us¬ 
ing the standard ISDC offline scientific analysis (OSA) software 
version 10.0. IBIS/ISGRI images for each ScW are generated in 
the 18-60 keV energy band. These images are used to produce the 
long-term light curve on the ScW timescale. The total spectrum is 
obtained using mosaic images as stated in the IBIS Analysis User 
Manutd- To extract the spectra we use the following parameters: 
size=l; posmode=-l; widthmode=-l. 

PSR B1259-63 was detected with a significance of 3.7 sigma 
in 3-35 keV combining JEMX1 & JEMX2 data. The low signifi¬ 
cance hindered further investigation of JEMX data. The source was 
detected at the 11 -a confidence level in 18-60 keV energy range 
by IBIS/ISGRI. The resulting 18 -60 keV light curve with a 3 day 
binsize is shown in Figure [5] The exposure time of each point is 
given in Table [6] The light curve shows a graduate decay of the 
flux with no striking variability. Please note that some of the flux 
detected from the position of PSR B1259-63 is actually coming 
from the nearby source 2RXP J130159.6-635806 , though as fol¬ 
lowed from the spectral analysis described below, the total 18-60 
keV flux of 2RXP J130159.6-635806 is less than a half of the flux 
of PSR B1259-63 in this energy range. 

The IBIS/ISGRI angular resolution (12 arcmin) does not al¬ 
low us to resolve PSR B1259-63 and 2RXP J130159.6-635806. 
Thus similar to the reconstruction of the hard spectrum of Suzaku 
one needs to fit the INTEGRAL spectrum as a sum of two sources. 
The quality of the INTEGRAL data allow us to fit the observed 
points with a single power law and the fitting procedure tends 
to attribute all the observed flux to PSR B1259-63. To prevent 


2 See http://www.isdc.unige.ch/integral/analysis for more information. 



Figure 5. PSR B1259-63 18-60 keV light curve with a 3 day binsize, as 
observed by IBIS/ISGRI. 


Table 6 . Integral observations ofPSR B1259-63. 


t t p 
(days) 

f(18 - 60 keV) 
(cts/sec) 

Exposure 

(ks) 

55 

0.90± 0.14 

27.21 

58 

0.81± 0.16 

23.65 

61 

0.67± 0.16 

23.27 

64 

0.57± 0.16 

25.42 

67 

0.50± 0.13 

40.27 


this we set a lower limit for the possible 20 - 60 keV flux from 
2RXP J130159.6-635806, to F minjaP = 10 11 erg s' 1 cm' 2 . This 
minimum value is about a half of the flux registered by NuSTAR in 
May 2015. With such a constraint we get a more realistic fit with 
X 2 = 9.93 for 9 degrees of freedom. The resulting 20 - 60 keV 
flux from 2RXP J130159.6-635806 varies from 10' 11 erg s' 1 cm' 2 
to 1.4 x 10' 11 erg s" 1 cm' 2 , and for the PSRB 1259-63 we get 
the following parameters: F 2 o -60 = 2.13/jjM x 10' 11 erg s' 1 cm' 2 , 
T = 1.67 ± 0.27. These parameters are in a good agreement with 
Swift observations taken during the same period (see Figure[7}. 
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3.5 The X-ray light curve 

Panel (b) of Fig. |T| shows the X-ray light curve of the system (note 
that the given 1-10 keV flux is not corrected for the absorp¬ 
tion). Observations made with different instruments at close orbital 
phases are consistent with each other, demonstrating good intercal¬ 
ibration. 

From data obtained at orbital phases similar to the archival 
observations of previous periastron passages, one can see that the 
system orbital light curve is stable over a time scale of several 
years. The stability of the orbital light curve allows us to use old 
and new data simultaneously while analysing the orbital evolution 
of the flux. Around the 2014 periastron passage our X-ray obser¬ 
vations started around the time of the first peak in the X-ray light 
curve at t ~ f p - 15 (Swl3,NStl). Similar to previous observa¬ 
tions we observed the gradual decrease of the flux, as the pulsar 
moved towards the periastron passage. Dense coverage with Swift 
after the periastron passage allowed us, for the first time, to observe 
the second rise of the flux in detail, and to determine that the sec¬ 
ond X-ray peak is almost twice as high as the first, similar to the 
behaviour commonly observed in radio. The general stability of the 
X-ray light curve allows us to use Chandra 2007 data to conclude 
that during the second rise the flux doubles in less than three days, 
and reaches the maximum less than in a week. After that, the X-ray 
flux gradually decreased in a good agreement with previous ob¬ 
servations. Flowever, detailed Swift observations show for the first 
time that start of GeV flare corresponds to the sharp change of the 
rate of the X-ray flux decrease, and also to a significant hardening 
of the X-ray spectrum, as discussed in Section [3~6l 


3.5.1 Fast variability in the light curve 


The fast, ~ 1 d, var iability of the system was previously suggested 
in lTam et alj J2015I ). This may indicate that the cooling of electrons 
in the PSR B1259-63 system happens already on daily time scales, 
which challenges the conventional interpretation of the decaying 
branch of the source light curve in terms of the gradual cooling 
of the injected particles. In order to investigate this possibility, we 
have searched for such fast variability in Swift and NuSTAR 2014 
data. 


The NuSTAR data alone, due to their sparse coverage 
of the source light curve, do not allow us to look for vari¬ 
ability on time scales longer than * 30 ksec ( the duration 
of the single NuSTAR observ ation). Ap plyin g the_ S truc 
ture Function ( SF) analysis 1 SimonettijCordes & Heeschei 


198j| 


Favre, Courvoisier & Paltanil 


t 


schen 

2005; 


Savchenko, Neronov & Courvoisietl 1201 2l : lYovk & Neronovl 


2013) to the light curves, we found no significant variability on the 


0.5-30 ksec time scale. 

Much better coverage of the periastron passage is given by 
the Swift observations, which provide light curve coverage in » 
10-30 ksec time bins, separated by ~ 1 d time intervals, which 
creates substantial di fficulties fo r the S F analysis. As a way to over¬ 
come this difficulty, iTam et aO d2015t) have fitted the light curve 
with the exponential decay model, demonstrating significant resid¬ 
uals from the fit. However, the intrinsic shape of the decaying 
branch of PSR B1259-63 light curve may deviate from the sim¬ 
ple exponential decay form, e.g. due to the movement of the emit¬ 
ting region through the complex environment of the system and the 
subsequent change in the cooling rate. 

In order to avoid complications due to the uncertain intrinsic 
light curve shape, we have chosen a different approach. In order to 


assess the variability at the time scale r, we first created an aver¬ 
aged light curve, using r as the averaging scale. This averaged light 
curve contains n av = (T max - T min )/r time bins. We then compute a 
difference of the original and the averaged light curves, which then 
can be quantified in terms of the^ 2 with n orig - n m . degrees of free¬ 
dom. Clearly, a very small value of r in this case would result in the 
averaged light curve perfectly matching the original light curve, re¬ 
sulting in the zero ^7 At the same time r = T max - T min would be 
equivalent to the fitting the entire light curve to the constant, which 
in the case of the PSR B1259-63 observations would result in a 
very large value of^ 2 , indicating significant variability. 

We repeated this procedure for the values of r in the range 
from 0.5 to 10 days with a step of 0.1 day. The result of this scan 
is shown in Fig. [6] This figure shows that the variability of the 
source starts to be significant above r as 4 days, corresponding to 
the rise/decay time scales of the main flares in the light curve. Thus, 
our analysis disfavours the daily variability of PSR B1259-63 sug¬ 
gesting that its light curve is rather smooth. 


3.6 Spectral Analysis 

The X-ray spectral analyses were done with NASA/GSFC XSPEC 
vl2.7.1 software package. A simple power law with a photoelectric 
absorption describes the data well, with no evidence for any line 
features. In Table [7] we present the results of the three parameter 
fit to the Swift and NuSTAR data. The uncertainties are given at 
the lcr statistical level and do not include systematic uncertainties. 
The graphical representation of the spectral parameters is shown 
in Figure [7] In Figure [7] new observations are shown along with 
the ones from previous periastron passages. The example of the 
NuSTAR spectra (NSt2 and NSt3) are given in Figure[3] 

The average hydrogen column density over the entire Swift 
2014 campaign is nn - 0.588 ± 0.027xl0 22 cnT 2 , estimated as 
a weighted mean of the values in Table [7] A similar analysis of 
the combined XMM-Newton and Suzaku 2007 data set suggests 
a lower value of n H = 0.487 ± 0.007xl0 22 cnV 2 . The two values 
are compatible with each other at the 3.6<x confidence level. When 
compared to the averaged value n H = 0.493 ± 0.007xl0 22 cnT 2 , the 
significance of the n H variation is reduced to 3.2<x. 

The dotted line in the Figure [7] shows the moment when the 
GeV emission was for the first time detected in the GeV band on 
a day time scale (31 days after the periastron passage, Caliandro et 
al. (2015, in press)). It is clearly seen that the rate of the source de¬ 
cay and its spectral index significantly changes after the flare. The 
decay becomes much more shallow, while the emission becomes 
much harder, with an average value of T = 1.5, instead of T = 2 
at periastron. No significant change of the column density has been 
observed. 


4 FERMI-FAT OBSERVATIONS AND RESULTS 

The Fermi- LAT results included in this paper are adopted from 
Caliandro et al. (2015, in press). The analysis of Fermi-LAH data 
was performed using the Fermi Science Toolfl 09-34-01 release. 
For the 2010 and 2014 periastron passages the analysis was car¬ 
ried out with Pass 7 reprocessed data belonging to the SOURCE 
event clasfl All gamma-ray photons within an energy range of 


3 http://fermi.gsfc.nasa.gov/ssc/ 

4 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Pass7REP_ 
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Figure 6. Left: false alarm probability of the variability in the PSR B 1259-63 light curve from the Swift data, see text for details. The 1,3 and 5 cr confidence 
levels are plotted with the dot-dashed, dashed and solid black lines correspondingly. Right: The light curve ofPSR B1259—63 over the 2014 observational 
period in the 1-10 keV band, shown together with its averaged version. The averaged light curve is produced with r = 4 days, which corresponds to the 
minimal scale, where the false alarm probability crosses the 3cr significance limit. 




Table 7. Spectral parameters for 2014 observations of PSR B1259-63. 


Data 

f-fp. 

F(l-10 keV) 

F^d-10 keV) 

r 

A h 

C* 

ndof 

Set 

d 

10 -11 erg cm -2 s -1 

10“ 11 erg cm -2 s -1 


10 22 cm - 2 



Swl3 

-14.56 

2 - 83 Td 

3 - 14 o n 

2-42 + “'!= 

1.63+0.21 

0.61+0.21 

189.95 

240 

Swl4 

-14.24 

7 ^+o:io 

2.07+0.16 

0.84+0.16 

260.99 

329 

NST1 

-14.19 

2.84 + ° : ® 

2.46 4 i 

1.84+0.02 

0.40+0.26 

0.90 

330 

Swl5 

-6.78 

2 - 14+ 006 

1 - 84+ ool 

1.71+0.12 

0.59+0.12 

376.91 

402 

Swl 6 

0.22 

L70+ O06 

1.32 + ||? 

2.04+0.15 

0.75+0.14 

295.15 

344 

NST2 

0.005 

i oo+8i8t 

1 77 +8:8^ 

1.95+0.03 

0.64+0.34 

0.98 

256 

Swl7 

2.16 

1.79« 


1.99+0.26 

0.57+0.24 

147.61 

190 

Swl 8 

2.35 

LM* 

1.30 + |jlf 

2 . 12 + 0.22 

0.97+0.22 

213.93 

247 

Swl9 

4.42 

] -37+OTO 

1.26 4 ! 

1.59+0.28 

0.20+0.25 

120.14 

155 

Sw20 

12.86 


2 - 40+ o I 9 

2.38+0.61 

1.19+0.61 

47.14 

63 

Sw21 

15.12 

5.26 + || 

3.66 4: 

2.16+0.19 

1.10+0.19 

237.50 

302 

Sw22 

17.70 


4.38 4 

1 . 86 + 0.10 

0 . 68 + 0.10 

457.72 

529 

Sw23 

20.76 

4.0lf{ 

•3 ^ 7 + 0.10 

1.85+0.11 

0.63+0.11 

403.67 

487 

Sw24 

23.96 

o C7+O.O9 

,,,78:18 
' 3 - -0.10 

1.57±0.11 

0.52+0.11 

413.16 

499 

NST3 

24.02 

q 77+0.03 

2.67 + “i 

1.61+0.01 

<0.67 


355 

Sw25 

29.19 

2 . 0 lf 

1.51+0.21 

0.39+0.22 

214.99 

313 

Sw26 

29.33 

2 .° 2 “ 

1.90 

1.48+0.26 

0.08+0.30 

220.49 

247 

NST4 

29.41 

2 . 86 ™ 

, ,,+8:Af 

2- 28 -0.09 

1.52+0.01 

<0.81 


336 

Sw27 

33.36 

77 ,;8$ 
z, ‘ / -0.09 
2.64 + ™ I 

1.72±0.14 

0.68+0.15 

338.63 

430 

Sw28 

37.56 

2.36 + SS| 

1.52+0.11 

0.48+0.11 

355.16 

496 

Sw29 

38.73 

2 42 + 8:?§ 
Z ^ Z -0.14 

2 07 + ^' 
z - w/ -0.16 

1.71±0.26 

0.54+0.24 

191.01 

250 

Sw3() 

40.11 

2.25 + °'n 

1 . 86 + “j| 

1.75+0.23 

0.67+0.22 

237.73 

319 

Sw31 

41.11 

1.84 + | j 

1.63-8: i 

1.62+0.28 

0.42+0.26 

201.54 

261 

NST5 

41.31 

2 ^j+ 8 :Al 

-0.Q2 

1 '9 8 —nips 

1.52+0.01 

<0.44 


343 

Sw32 

41.55 

7 7 ^+ 0 .17 

2 . 22 + jj'?„ 

1.09+0.28 

0.15+0.29 

127.01 

203 

Sw33 

42.85 

2,4 

, q ,+8:i8 

1.68+0.18 

0.50+0.17 

281.04 

354 

Sw34 

44.52 

2. 1 8 + q' j j 

2-f> 4+ |'l 

1.41+0.25 

0.26+0.23 

180.19 

248 

Sw35 

51.60 

7 07+O.O6 

-0.Q7 

2 io +01( J 
-Q-iQ 
1.96 + 9'9Z 

7 17+8:45 
-0.08 

1.23+0.13 

0.49+0.14 

371.76 

418 

Sw36 

Sw37 

52.78 

53.68 

1 76 +011 
-QQ9 
1.73 + °® 

1.63+0.21 

1.40+0.16 

0.72+0.22 

0.61+0.18 

210.88 

267.29 

244 

335 

Sw38 

58.71 

L96+ ooi 

L59+ 006 

1.72+0.15 

0.83+0.16 

327.19 

358 

Sw39 

63.02 

L56+ 00S 


1.62+0.16 

0.76+0.17 

262.48 

338 

Sw40 

63.52 

! 41 78:8? 

1>41 -0.05 

, ,+8:88 
1 ’ Ji -0.05 

1.26+0.16 

0.40+0.17 

256.06 

325 

Sw41 

64.85 

1 OQ+0.05 

_-004_ 

1 1 n+0.05 

_ 1 -0.05_ 

1.53+0.16 

0.72+0.17 

271.68 

343 


* Relatively poor statistics in some of the Swift spectra forces us to use C statistics in data analysis. 
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Figure 7. PSR B1259-63 orbital evolution of (1-10 keV) light curve (top panel), spectral index (middle panel) and hydrogen column density (bottom panel), 
as seen with Swift and NuSTAR during the 2014 periastron passage along with the old XMM-Newton Suzaku and Swi/lobservations. 1-10 keV flux of the 
source is given in units of 10 11 erg cnT 2 s 1 . 


0.1-100 GeV and within a circular region of interest (ROI) of 10° 
radius centered on PSR B1259-63 were used. To reject gamma- 
ray contamination originating from the Earth’s limb, we selected 
events with zenith angle < 100°. The gamma-ray flux and spectral 
results of PSR B1259-63 presented in this work were calculated 
by performing a binned maximum likelihood fit using the Science 
Tool gtlike. The spectral-spatial model constructed to perform the 
likelihood analysis includes Galactic and isotropic diffuse emission 
components as well as known gamma-ray sources within 15° of 
PSR B1259-63 based on a catalogue inte rnal t o the Fermi collab¬ 
oration ( (now released as 3FGL. lAcero et all 120151) ). The spec¬ 
tral parameters were fixed to the catalogue values, except for the 
sources within 3° of PSR B1259- 6fl For these latter sources, the 


5 http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/source. 


flux normalization was left free. PSR B1259-63 itself was mod¬ 
elled as a single power-law with all spectral parameters allowed to 
vary. The contribution of the Galactic and isotropic diffuse emis¬ 
sions within the analyzed ROI was estimated performing a prelim¬ 
inary maximum likelihood fit. For the 2010 periastron event, this 
fit included data from 5 months prior to 1 year after periastron. For 
2014, this fit included data from 5 months prior to periastron. The 
resulting scale factors of the isotropic and Galactic diffuse tem¬ 
plates were kept fixed in the analysis. 

The 2014 Fermi -LAT observations showed a persistent nature 
of the GeV flare with the onset of the flaring activity starting 31 
days after periastron (Caliandro et al. (2015, in press)). It turned 
out that the 2014 GeV flare exhibits a similar average flux level 
and spectral shape with the flare of 2010. The details of the 2014 
GeV flux evolution, differs however from the 2010 light curve, see 
_podell6a)itifi]Figure [T] The 2010 flare has a higher peak flux, and a 
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rapid decreasing evolution (Figure 2a). The 2014 flare has a peak 
flux lower by a factor of 1.6, which then persists rather than falling 
rapidly. These differences may be due to inhomogeneities in the 
shape, density, or extent in the circumstellar disk of the Be star 
(e.g. the equivalent width of the Ho- line before the periastron was 
somewhat higher in 2010, than in 2014, see panel (c) of Figure Q} 


5 DISCUSSION 


The broad band non-thermal emission from PSR B1259-63 is pro¬ 
duced by high-energy particles ac celerated at the shock formed at 
the interface of the pulsar winds dTavani & Arons! fl997t) . The or¬ 
bital modulation of different components of the broad band spec¬ 
trum is generally understood as being due to the orbital variabil¬ 
ity of geometrical parameters, such as the geometry of the contact 
surface of the stellar and pulsar winds. The overall "two-bump" 
variability pattern of the source, evident from Fig. [7] is generally 
explained by the passage of the pulsar through the dense equatorial 
wind of the Be star. This happens twice per orbit, just before and 
just after the periastron. 

The most puzzling feature of the orbital modulation of the 
source is t he GeV band flaring activity discovered by Fermi/LAT 
telescope dAbdo et aPl201 ll : |chernvakova et alj|2014l) . The flare 
occurred unexpectedly a month after periastron and ten days af¬ 
ter the post-periastron passage of the pulsar through the equatorial 
wind of the Be star. By this time, the pulsar is already exiting from 
the dense equatorial stellar wind of the Be star and the overall ac¬ 
tivity of the source is decreasing in all energy bands, from radio to 
X-ray. The sudden brightening of the source at this moment of time 
was not expected before the Fermi observations. 

The GeV flaring appeared even more puzzling in the absence 
of clear counterparts to the flare in other energy bands. The flare 
appeared as an "orphan" increase of the GeV flux, not accompanied 
by changes in the source state in other wavebands. 

The origin of this "orphan" flare was w i dely discussed in 
the lit e rature (see e.g. IP etri & Dubusl d201ll): iKha ngulvan et alj 


d2012l) : iMochol & Kirkl d2013l) : iDubus & Ceruttil d2013l) ). How- 


ever, no consensus on the nature of the flare was reached, because 
of the significant lack of information on the flare properties. In par¬ 
ticular, before the new observational campaign around the 2014 pe¬ 
riastron passage, it was not known if the GeV flare was a unique 
event or if it is systematically repeating at each periastron passage. 
Even if the flare is recurrent from orbit to orbit, it was not clear 
if it always occurs at the same orbital phase or if it could occur at 
different phases. It was also unclear if there is only one flare per pe¬ 
riastron passage or if the flaring activity could appear as a sequence 
of flares at different orbital phases. 

New observations in the GeV band in 2014, combined with the 
multi-wavelength observations reported here provide the missing 
observational constraints necessary to clarify the nature of the flare. 
It is now clear that the GeV flaring is recurrent from orbit to orbit. 
It always occurs at one and the same orbital phase, some 30 days 
after the periastron passage. 

X-ray observations provide a detailed picture of the behaviour 
of the source throughout the GeV flaring period and show that the 
GeV flare is not "orphan". It is accompanied by a specific hard state 
of the source in the X-ray band (see discussion in Section [T6l i. 

Optical spectroscopy data reported here also demonstrate that 
the moment of the onset of the flare is not random. It coincides with 
a sudden decrease of the equivalent width of Her line, which char¬ 


acterises the state of the equatorial disk of the Be star (see Section[2] 
and panel (c) of Fig. |TJ- 

These new ob servations confirm a mode l of the flare which 
was put forward bv lchernvakova et alj J2014I ). Within this model, 
the GeV flare occurs at the moment of partial destruction of the 
equatorial disk of the Be star by the passage of the pulsar. Estimat¬ 
ing t he mass of the Be star di sk from the equivalent width of Her 
line JChernvakova et al1l2014l) . one could find that at the moment 
of the GeV flare, the disk mass has decreased by a factor of five, 
from about 2 x 1CT 8 M 0 down to 4 x 10~ 9 M G , within just five days 
after the onset of the flare. The nature of such a dramatic event of 
destruction of the Be star disk has to be further investigated, but the 
consequences of the event for the properties of non-thermal emis¬ 
sion from the system could be understood in a straightforward way. 

Destruction of the equatorial disk also destroys the well- 
organized geometry of the interacting pulsar / stellar winds sys¬ 
tem. The regular bow -shaped contact surface of the two winds 
dTavani & Arons! 1 19971) is destroyed. Instead, a chaotic system of 
clumps originating from the dense Be star disk produces a contact 
surface of highly irregular geometry. 

The destruction of the regular contact surface also closes the 
escape path for the unshocked pulsar wind which was previously 
able to escape along a cone which points away from the Be star. The 
pulsar wind which could not escape releases all its power inside the 
system. This leads to the increase of the luminosity of the source, 
up to the ~ 68% of the _spin-down power of the pulsar in 2010 
and ~ 50% in 2014 l lAbdo et all d201 ll) . Caliandro et al. (2015, in 
press)). 

Fig. [8] shows a comparison of the new measurements of the 
spectrum of the flare and pre-flare states of the source with the 
model of lAbdo et alj d201 ll) . The Fermi/LAT data from the both 
observed flares are combined, which allows for a measurement of 
th e y-ray spec trum down to lower and at higher energies, compared 
to lAbdo et al .1 d201 ll ). Higher quality NuStcir data reveal the X-ray 
counterpart of the flare. 

These higher quality data are still consistent with the model in 
which the flare emission is attributed to the synchrotron emission 
from electrons injected with a close to monoenergetic spectrum and 
cooled by the synchrotron energy loss. The synchrotron cooling 
leads to the formation of a dNJdE oc E~ 2 low energy tail of the 
electron spectrum. Synchrotron emission from electrons in the tail 
has the hard spectrum dN y /dE oc E^ 15 observed in the X-ray band 
throughout the flare. 

Neither the energy nor the origin of the nearly monoenergetic 
electrons could be firmly established based on the observational 
data. This is because the maximal energy of the synchrotron emis¬ 
sion is close to the self-regulated cut-off in the synchrotron spec¬ 
trum at E cu , ^ m e c 2 /a ^ 100 MeV, where a is the fine structure 
constant. Such a cut-off occurs when electrons are accelerated at 
the maximal rate dE/dt ^ eB and balanced by the synchrotron loss 
rate. In this regime the maximal energy of the synchrotron emission 
does not depend on either the magnetic field B or on the energy of 
electrons. Assuming that the magnetic field in the pul sar a nd stellar 
wind interaction zone is about B ~ 1 G, iTavani& Arons! l l 19971) 
provides an estimate of electron energy E e ~ 10 14 eV, which is 
necessary to produce synchrotron emission in the 100 MeV energy 
range. 

Contrary to the flaring component, the persistent component 
of the broad band spectrum does not have a sharp peak at 100 MeV. 
Instead, the X-ray-to-GeV spectrum is flat in the EF e representa¬ 
tion, i.e. it has the form dN/dE oc E~ 2 . Synchrotron emission with 
such a spectrum is produced by electrons with a powerlaw spectrum 
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Figure 8. Broad band spectrum of the source for the flaring state (black 
data points) and for the period of periastron passage excludin g the GeV 
flare (b lue data points) superimposed on two model spectra from, Abdo et ah 
1201 ill for the synchrotron (solid lines) and inverse Compton (dashed and 
dotted lines) emission during the GeV flare (grey) and during the rest of the 
periastron passage (light blue) periods. 


dNJdE oc E~ 3 . The synchrotron spectrum hardens below keV and 
the slope of the radio-to-X-ray spectrum is dN/dE oc E~ 15 , which 
corresponds an electron spectrum with a slope dN e /dE oc E~ 2 . 
Thus, the electron spectrum responsible for the persistent emis¬ 
sion throughout the two passages of the Be star disk is a broken 
powerlaw (dNJdE oc E~ r ) with the slope softening from T =* 2 
to T — 3 at the energy E e - 1 [6/1 G]^ 1/2 TeV. The softening by 
AT — 1 is typical for the effect of the synchrotron cooling on the 
spectrum. The synchrotron cooling modifies the electron spectrum 
in the energy range above 1 TeV and it fails to do this below 1 TeV. 
The absence of cooling in the energy band below 1 TeV could be 
attributed to the escape of the sub-TeV electrons from the system. 
The synchrotron cooling time is t s - 400 [6/1 Gp 2 [E e /1 TeVp 1 s. 
It is comparable to the escape time t esc = R/s - 300 [6/10 13 cm] s 
for a region which is approximately the size of the binary separa¬ 
tion distance R ~ 10 13 cm. 

The injection spectrum of electrons responsible for the persis¬ 
tent emission during the disk passages, dNJdE oc E~ 2 is differ¬ 
ent from the spectrum of electrons injected during the GeV flare 
(which is close to mono energetic). This suggests that the two pop¬ 
ulations of high-energy electrons are produced by different acceler¬ 
ation processes and/or originate from different sources. The natural 
possibility for different emission sites is electrons originating from 
the pulsar wind and electrons from the stellar wind. The two pop¬ 
ulations of electrons might also appear because of the presence of 
two different acceleration sites in the system: the non-relativistic 
shock from the stellar wind side of the contact surface and rela¬ 
tivistic shock at the pulsar wind side of the winds contact surface. 

In an alternative model of the flaring activity of the source 
l iAbdo et al.l201 ll) . the "orphan" flare in the GeV band was assumed 
to be produced via the Bremsstrahlung mechanism. A flare would 
correspond to a sudden increase of the density of the medium, 


which could be the debris of the destroyed equatorial disk of the 
Be star. Such a model is disfavoured by the observations of the X- 
ray counterpart of the flare. The Bremsstrahlung emission produces 
a negligible contribution to the X-ray flux and could not be respon¬ 
sible for the hard spectrum flaring activity coinciding with the GeV 
flare period. 

Within the synchrotron scenario, the flaring activity of PSR 
B1259-63 is similar to the flaring activity of the Crab pulsar where 
the GeV band flares are due to the synchrotron emission from the 
highest energy electrons accelerated at maximal possible rate. A 
study of the comparison of the detailed properties of the GeV flares 
in the two sources could potentially be useful for understanding the 
origin of the high-energy electrons and the location of the high- 
efficient acceleration sites inside the two sources. Observations re¬ 
ported in this paper show that the triggering of the flares in PSR 
B1259-63 occurs simultaneously with an instability which destroys 
the equatorial disk of the Be star companion to the pulsar. It is not 
clear if it is the disk instability which triggers an instability of the 
pulsar wind and of the interface of the pulsar and stellar winds, or 
the causal link is in the opposite direction. In the case of Crab, the 
mechanism of triggering the flares could only be an intrinsic insta¬ 
bility of the pulsar wind or of the interface of the pulsar wind and 
external medium. An alternative possibility is not available because 
of the absence of the companion massive star. It remains to be seen 
if the flares of Crab and PSR B1259-63 are triggered by the same 
type of instability and are powered by the same type of acceleration 
mechanism. 
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